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Scattering of shock waves in QCD 
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The cross section of heavy-ion collisions is represented as a double functional integral with the sad- 
dle point being the classical solution of the Yang-Mills equations with boundary conditions/sources 
in the form of two shock waves corresponding to the two colliding ions. I develop the expansion 
of this classical solution in powers of the commutator of the Wilson lines describing the colliding 
particles and calculate the first two terms of the expansion. 
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I. INTRODUCTION 

Viewed from the center of mass frame, a typical high- 
energy hadron-hadron scattering looks like a coUision of 
two shock waves (see Fig. 1). Indeed, due to the Lorentz 
contraction the two hadrons shrink into thin "pancakes" 
which coUide producing the final state particles. The 
main question is the field/particles produced by the col- 
lision of two shock waves. On the theoretical side, this 
question is related to the problem of high-energy effective 
action and to the ultimate problem of the small- a; physics 
- unitarization of the BFKL pomeron and the Froissart 
bound in QCD (see Q, 0). On more practical terms, 
the immediate result of the scattering of the two shock 
waves gives the initial conditions for the formation of a 
quark-gluon plasma studied in the heavy-ion collisions at 
RHIC (see e.g. the review Q)- 

The collision of QCD shock waves can be treated us- 
ing semiclassical methods. The basic idea is that at high 
energy the density of partons in the transverse plane be- 
comes sufficiently large to give the hard scale necessary 
for the application of perturbation theory 0, Q . The ar- 
guments in favor of this are based on the idea of parton 
saturation at high energies H, 0- Consider a single 
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FIG. 1: High-energy scattering as a collision of two shock 
waves. 
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shock wave - hadron, moving at a high speed (in the cm. 
frame). The energetic hadron emits more and more glu- 
ons and the gluon parton density increases rapidly with 
energy. This cannot go forever - at some point the recom- 
bination of partons balances the emission and partons 
reach the state of saturation with the charactristic trans- 
verse momenta (the "saturation scale") being Qs ~ e'^^ 
where 77 is the rapidity H [H [H . Such an energetic 
shock wave with large density of color charge is called the 
Color Glass Condensate 0, . 

Within the semiclassical approach, the problem of sca- 
tering of two shock waves can be reduced to the solution 
of classical YM e quat ions with sources being the shock 
waves (see also |l4j). At present, these equations have 
not been solved. There are two approaches discussed in 
current literature: numerical simulations |l5j | and expan- 
sion in the strength of one of the shock waves 

Note that the collision of a weak and a strong shock 
waves corresponds to the deep inelastic scattering from a 
nucleus (while the scattering of two strong shock waves 
describes a nucleus- nucleus collision). In the present pa- 
per I formulate the problem of scattering of shock waves, 
find the boundary conditions for the double functional 
integral for the cross section and decribe the expansion 
in the commutators of two shock waves equivalent to the 
expansion in strength of one of the waves. The main 
technical result is the calculation of the second term of 
this expansion (the first term can be restored from the 
current literature) 

The paper is organized as follows. Sec. 2 is a more for- 
mal introduction: I outline the idea of the factorization 
of the hadron-hadron cross section into the formation of 
two shock waves and their scattering. In Sec. 3 I discuss 
the rapidity factorization and define what is a scattering 
of QCD shock waves. In Sec. 4 I find the Lipatov vertex 
of the gluon emission and in Sec. 5 reproduce the kx- 
factorization valid in the first order in the commutator 
expansion (for the pA scattering). In Sec. 6 I obtain the 
first-order effective action and reproduce the non-linear 
equation for the small-x evolution of Wilson lines. Sec. 
7 outlines the calculation of the second-order classical 
field in the while the details of the calculation are given 
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FIG. 2: Integration over the central rapidity gluons in a semiclassical approximation. 



in the Appendices A-C. The exphcit form of the vertex 
of gluon emission by two Wilson hnes in the shock-wave 
background is presented in the Appendix D. 



II. BASIC IDEA: TWO-STEP INTEGRATION 
OVER RAPIDITY 



In this section I outhne how the hadron-hadron colK- 
sion at high energy is related to the scattering of shock 
waves. The basic idea of the approach is the two-step in- 
tegration over rapidity in the double functional integral 
for the cross section. At first, let us define this integral. 



A. Double functional integral for the cross section. 

A total cross section is a product of an amplitude and a 
complex conjugate amplitude so the functional intergral 
for the cross section has double set of fields: to the right 
of the cut and to the left of the cut. A typical functional 
integral has the form: 



DADAJ{pA)J{pB)e-'^^^'^ J{-Pa)J{-Pb) e^^'^^ 

(1) 

where the currents J {pa) and J (job) describe the two 
colliding particles (say, photons). Throughout the pa- 
per, fields to the left of the cut will be represented by the 
calligraphic letters while those to the right of the cut by 
usual letters. The boundary conditions are such that the 
fields A and A coincide at i — s- cx), reflecting the sum- 
mation over the final states implied in the definition of a 
total cross section. The propagators for such functional 



integral reproduce the Cutkovsky rules (cf. ref. [l^) : 



ah 



{Al{x)Al{y)) ^ g,J^ 
{Al{x)Al{y)) ^ g,J^^' 



{Al{x)Al{y)) ^ -9,J'^' 



167r4 



-ik{x-y) _ 



167r4 



i 



(2) 



167r4 



fc^ — ie 

,-ifc(x-j/)r,_X/';,2 



We are interested in the number of gluons produced 
per unit rapidity which is given by the average of the 
creation operator over the final state {a{k±,rj)a^ {k±,r])) 
(see the discussion in ref. [23|). In terms of functional 
integrals this can be rewritten as 

ng{k^,Tj) = lim [dADAJ{pa)J{pb) (3) 
X e-^^f-^) k^Ank)k^A^{-k)e'^^^h{-pA)J{-PB) 

Throughout the paper, the sum over the Latin indices 
runs over the two transverse components (while the 
sum over Greek indices runs over the four components 
as usual). As we shall see below, the Lipatov vertex 
of gluon emission i?^(/c) = limfc2^g ^^(/c) is transverse: 
fc^i?^(/c) —) so we can replace AiAi in Eq. Q by the 
sum over all four indices 



n, 



g{k^,Tl) = - lim I DADAJ{pa)J{pb) (4) 
X e~'^(-^) eAl{k)k' A^'{-k)e'^^^^ J{-pa)J{~Pb) 



B. Two-step integration. 

The integration over the gluon fields in the functional 
integral (QJ will be done in two steps according to the 
rapidity of the gluons. Let us introduce two rapidities 
•qi and r/2 such that i]a > t]i > V2 > Vb- Consider a 
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typical Feynman diagram for the gluon production Q 
shown in Fig. |21 At first, we integrate over the fields in 
the central range of rapidity 771 > > 772 and leave the 
fields with 77 > 771 and 77 < 772 in the form of external 
shock waves. In the semiclassical approximation there is 
only one gluon emission described by the Lipatov vertex - 
Fourirer transform of the classical field at the mass shell. 
The result of the integration is the product of two Lipatov 
vertices which depend on that "external" fields. It is 
easy to see that Lipatov vertices depend on these gluon 
fields through Wilson lines - infinite gauge links ordered 
along the straight lines coUinear to 771 and 772. Indeed, 
in the target frame gluons with rapidities 771 > 77 > 772 
are very fast so their propagators in the background of 
"target" gluons reduce to the gauge factor ordered along 
the straght line classical trajectory. 

Thus, in the semiclassical approximation we get (see 
Fig. H 

ng{kr:V) = j DADADBDBJAJAe''^^^^+'^^'^^ (5) 
where R{k]V,U ,V ,U) \s a. product of two Lipatov vertices 



R{k;V,U,V,U) 



Rf,{k;V,U)n''{-k,V,U) (6) 



The Lipatov vertex R^{k; V, U) = hmfe2_^o k'^A^{k; U, V) 
is an amplitude of the emission of a gluon with momen- 
tum k by the two Wilson lines U and V. It depends on 
the gauge, but the product of two Lipatov vertices © is 
gauge-invariant due to the property k^R^^lk) ~ 0. Here 
A, A are fields with rapidities 77 > 771 and S, B with ra- 
pidities ?7 < 772. The Wilson lines V and U are made 
from A and B fields, respectively: 



V{xj_) = [OOTT-I + X±, — 007ll + X±], 
U — [007l2 + X±, —00712 + X±] 



where ni and rt2 are the unit vectors corresponding to 
rapidities 771 and 772 while [x, y] is a shorthand notation 
for the straight-line ordered gauge linkconnecting points 
X and y: 



[x,y] = Pe^sIodu{x-y)''A^{ux+{l~u)y) 



(8) 



It may seem that the result (O depends on the ar- 
tificial "rapidity divides" 771 and 772. This dependence 
should be canceled by the gluon ladder on the top of 
the Lipatov vertices which is outside the semiclassical 
approximation. The common belief is that the all the 
evolution can be attributed to either upper or lower sec- 
tor: to calculate ng{k±, 77), we choose 771 and 772 such that 
both q;s(77i — 77), (77 — 772) ^ 1 so the part of the evolu- 
tion between 771 and 772 can be neglected and we can use 
the semiclassical aproximation for the functional integral 
(O over the region of rapidity 771 > 77 > 772 (see e.g. the 
review Q). Eventually, after solving the classical Yang- 
Mills equations for this functional integral, we can put 



771 = 77 in 1/ and 772 = 77 in U. The independent evolution 
in the upper (or lower) sectors leads to the parton satu- 
ration so the shock waves U and V have a form of Color 
Glass Condensate 

To find the Lipatov vertices one needs to solve the 
classical YM equation with the sources proportional to 
shock waves U and V . As we mentioned, these equations 
have not been solved yet and there are two approaches 
discussed in the literature: numerical simulations and ex- 
pansion in the strength of one of the shock waves. In this 
paper I develop the second approach in a "symmetric" 
way as an expansion in commutators [U, V] , and calcu- 
late the second term of the expansion (the first one can 
be restored from the literature). 



III. RAPIDITY FACTORIZATION AND 
SCATTERING OF THE SHOCK WAVES 

A. Rapidity factorization 

In this section we define the scattering of the shock 
waves using the rapidity factorization developed in 0, 
I2H . Consider a functional integral for the cross section 
Q and take some "rapidity divide" 771 such that 77^ > 
Vi > Vb- 

Throughout the paper, we use Sudakov variables 



k = api + /3p2 + k± 



(9) 



and the notations 



X, = p^Xf, = x+ = -^{x° + x^) (10) 



V2 
1 

72 



ix° - x^) 



(7) Here pi and p2 are the light-like vectors close to pA and 



Pb- Pa = Pi + ^P2, Pb ^ P2 

Let us integrate first over the fields with the rapidity 
77 > 771. From the viewpont of such particles, the gluons 
with 77 < 771 shrink to a shock wave so the result of the in- 
tegration is presented by Feynman diagrams in the shock- 
wave background. In the covariant gauge, the shock-wave 
has the only non-vanishing component A, which is con- 
centrated near = 0. A typical Green function G(x, y) 
at a;,, < in the background- Feynman gauge has the 
form [221 

{A{x)A{y)) = 

2 — 1 r 4i 

dzS{-z^){x\^ -\ 2ag^MzUz H {dJUzUz)p2u 

s p'^ — le ^ s 

_^^,)_ ^E^dliuU)Mz\^\y) (11) 
as^ J p'^ + 1€ 

where lAz = [oopi + z^, — oopi + z\_] is made form the left 
fields A while — [oopi -I- z^, — oopi + z^_] is made from 
As. 



'-fPi- 
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Similarly to the case of the usual functional integral 
for the amplitude, in order to write down factorization 
we need to rewrite the shock wave in the temporal gauge 
Aq = 0. In such gauge the shock-wave background has 
the form 



where 



9 



9 



(12) 



(13) 



are the pure gauge fields (fiUing the half-space < 0). 
Note that the choice (1121) is different from the choice Ai = 
Ui9{x.j,) adopted in [l4ll2^. The reason is the following: 
when we calculate the amplitude, it is natural to use 
the redundant gauge rotation to get rid of the field at 
t = —oo (althougth the choice t = oo \s equally possible). 
On the contrary, for a double functional integral of the 
form we have two independent integrations aX t = 
—oo and it is impossible to eliminate both of them. We 
can however gauge away fields at t = oo because of the 
boundary condition A — A\t=oc- (Stricltly speaking, we 
cannot gauge away all the fields; what we can do is to 
forbid a pure gauge fields t = oo which is enough for our 
purposes since it puts forward the choice p2|l over the 
choice Ui9{xf)). 

The Green functions in the background H12|) differ from 
those of (|ll|l by a simple gauge rotation. Their explicit 
form is presented in the Appendix I A II 

The generating functional for the Green functions in 
the Eq. H12|l background is obtained by the generaliza- 
tion of the generating functional of |22i| to the case of a 
double functional integral: 

DADAJ{pA)J{-PA)e-'^'-^'^+'^'^''^''°°-^'-°°^'"" 

^iS{A)-if cPz^(oo,F,,,ca)lU'^^ ^^4) 

where the Wilson-line operator 

(C50,.F^„(X3)'^ = (15) 

/oo 
du [oop2,Up2]z^t^i{up2 + Z±)[up2,-OOp2]z, 
'OO 

is made from the "left" fields and the operator 

{oo,F^,,ooY = 

2tr r / du [oop2,up2]zF^i{up2 + z±)[up2, -oop2]z 



from the "right" fields A^. 

It is easy to see that the functional integral (|14|l gener- 
ates Green functions in the Eq. (|12|) background. Indeed, 
let us choose the gauge A^ ~ A^. — for simplicity. In 
this gauge (oo, .f",^, oo)" = ^i(oop2 + z^) - Ai{-oop2 + 



zj_) and (00,^,^,00)" = Ai{oop2 + z±) - Ai{-oop2 + zj 
so the functional integral H14I) takes the form 



J DAD AJ (pa) J {-pa) 



(16) 



X e 



X e 



-iS(A)+iJ d^z±(Ai{oop2+z±}-Ai{-oop2+z±)yu" 
iS{A)-iJ d^z±{Ai{oc3p2+z±)-Ai{-cx>p2+z±))'^U" 



Let us now shift the fields Ai Ai+Ai and Ai Ai+Ai 
where Ai = Ui9{—x^) and Ai = lAid{—x^,). The only non- 
zero components of the classical field strength in our case 
are F,i = —Ui6{-x^,) so we get (for the "right" sector) 



S{A + A) 



X A'Fj 



- J d^zD'F^.A* - - J d^zj_dz, 

-f iA^(5v - 2»^M^M'' + - (17) 



In the A* = gauge the first term in the r.h.s. of Eq. 
(|17|l vanishes while the second term cancels with the cor- 
responding contribution ~ —(Ai(oop2 + z±) — Ai(—oop2 + 
z±))'^U°'^ coming from the source in Eq. H14|l . Similar 
cancellation occurs in the left sector so we get 



j DADAJ{pa)J{-Pa) 



(18) 



X e 



X e 



~iS{A)+iS d^z±_ [Ai{cop2+z±_)-Ai{-cop2+z±_)YUf 
iS(A)-iJ d'^z±[Ai(oap2+z±)-Ai{-oop2+z±)]'^U^^ 



J DADAJ{pa)J{-pa) 



X e 



-if zAf" [V^ g^i^-2if,„,)A'' zA" (D''- g,,u-2iF,„,)A'' 



which gives the Green functions in the Eq. H12|l back- 
ground. 

To complete the factorization formula one needs to in- 
tegrate over the remaining fields with rapidities 77 < ?/i : 

j DADAJ{pA)J{pB)e-'''^^^ e'^^''\ji-pA)J{-PB) 
= J DADAJ{pa)J{-Pa) j DBDBJ{pb)J{-Pb) 



X e 



-iS{A)—iS{B)+i^ d^z± (oc,JF^i ,cxd)^ (oo,C/^.i,oo)^e^e^ 



^ giS(A)+iS(B)-i/ 2:^(00, F^i,oo)"(cx;,G„i,cx;)"eje5' ^-|^g^ 

As discussed in [mill the slope of Wilson lines 
is determined by the "rapidity divide" vector e,,^ = pi + 
e~^^P2- (From the wiewpoint of A fields, the slope ei 
can be replaced by p2 with power accuracy so we recover 
the generating functional H14|) with Ui = (00, G*;, 00), 
= {00, Q^i, 00)). 

Applying the factorization formula H19|l two times one 
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gets: 



- lim / DADAJ{pa)J{pb) (20) 

= - J DADAJ{pA)J{^PA)e-""^-^~>+'^^^^ 
j DBDBJ{pB)J{~PB)e- 



-iS(B)+iS(B) 



X I DCDC lira p^C''^'{p)p^Cl{~p)el^Y> - iS{C) 

- [oo,A^j,cx)]"[oo,Cj.i,cx)]°^ +e^e2(^(cx),C^i,oo)" 
X (cx3,;Bi,i,oo)^ - (oo,C^i, 00)^(00, S,,^, 00)°^ I 

where the slope is ei = pi + e^^^p2 for the [...] Wilson 
lines and 62 = pi + e~^^p2 for the (...) ones. 

The functional integral over the central range of rapid- 
ity ?7i > 77 > 772 is determined by scattering of shock two 
shock waves: 

DCDC \im^k^C'"'ik)k^C^{~k)exp iS{C) - iS'(C) 

Jd'z^[ (yr [00 , 00]: - vr a, , 00]:) 

+ (z^r(oo,C^„cx3)^~C/f(oo,C^„oo):)e^}] (21) 

where = [00, ^^j, 00], = [00 , , 00] , = 
62(00, B^i, 00), and J7i = 62(00, 00) are the pure- 
gauge "external" fields (to be integrated over later). 
With a power accuracy 0(m^/s), we can replace ei by 
pi and 62 by p2 : 



£iCi:iCJjm^A:2c''^(fc)fc2c;^(-fc) exp {^^(C) - iS'(C) 

+ ijSz^ [00, C.„ 00] J + (00, C„, oo)^Zif 

- K"[oo,C.„oo]^- (oo,C„,oo)^t/r)} (22) 

The saddle point of the functional integral H22|l is deter- 
mined by the classical equations 



5C?. 



-{S{C) - tj d^z^(vr[^,C.,,^]t 
+ i^,C„,^)tUT)} = 
^{siC) - I j dV(vr[oo,C.„oo]? 

+ (oo,C„,oo):Z^r 



(23) 

At present it is not known how to solve this equations 
(for the numerical solution see ^3)- next section 
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FIG. 3; Particle production due to the scattering of the two 
shock waves. 



we will develop a "perturbation theory" in powers of the 
parameter [U,V] ^ g^[Ui,Vj] ^ g^[Ui,Vj]. Note that the 
conventional perturbation theory corresponds to the case 
when Ui,Vi ~ 1 while the semiclassical QCD is relevant 
when the fields are large {Ui and/or Vi ~ i). 



B. Expansion in commutators of Wilson lines 

The particle production due to scattering of the two 
shock waves in QCD (see Fig. O is determined by the 
functional integral H22I) (hereafter we switch back to the 
usual notation A^ for the integration variable and 
for the field strength) 

R{p; U, VM, V) = - ^hm^ j DAD A k^A'^^'ik) 

X k^A^^i-k) exp I - iS{A) + i J d^z^^Vf [00, J^,„oo]^ 

+ (00, oo)^Z^f ) + iSiA) - i J (fz^ (y^' 

X [00, F„, 00]^ + (00, F„, oo)^C/f } (24) 

Taken separately, the source Ui creates a shock wave 
Ui9{—x^) and Vi creates Vi9{—x,) (to the left of the cut, 
Ui generates the classical field Ui9{—x^) and Vt generates 
Vi6{—x,)). In QED, the two sources Ut and Vi do not 
interact so the sum of the two shock waves 



A 



(0) 



-.Ure{-x,) + v^e{-x,), A': 



(0) 



Af^ ^u,e{^x^) + ve{-x,), 



A 



(0) 



(25) 



is a classical solution to the set of equations H23|l . In 
QCD, the interaction between these two sources is de- 
scribed by the commutator g\Ui,Vk\ (the coupling con- 
stant g corresponds to the three-gluon vertex). It is nat- 
ural to take the trial configuration in the form of a sum 
of the two shock waves and expand the "deviation" of 
the full QCD solution from the QED-type ansatz H25|) 
in powers of commutators [[/, V]. To carry this out, one 



shifts A 



A + Af^ , A ^ A-^- AT' in the functional 



(0) 
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integral (|24() and obtains 



A. Pure gauge fields at t = — oo 



R(k;U,V,U,V) ^ - [dADA \im k^A''''(k) 

J fe^^O 

A^(-fc) exp {i jd^zi^ - {\a'"D^,A'' 
A'^D^^A- - r^A^ + r^A^) } , (26) 



X fc2 
1 



Here D^,^ = T)'^{A)g^y - 2iF^^ (P^^ = D'^{A)gf,u - 
2j^/ii/) is the inverse propagator in the background- 
Feynman gauge [s^ and I)iC?^ is the hnear term for the 
trial configuration p5|l 



= -D''F^ke{-z,)e{-z,) 



(27) 



+ ig[Ui, V] [pi^,9{z,)d{z,) - P2t.e{z,)5{z^) 

- -D^T^ke{~z^,)e{-z,) 

+ ig[U„V'](pi^e{z^)6{z,) - p2^0{z,)6{z,) 



where Fik = -ig[Ui, Vfe] - i ^ fc (and = ~ig[Ui, Vk] - 
i k). 

Expanding in powers of T{T) in the functional integral 
(|26|l one gets the set of Feynman diagrams in the external 
fields H25|) with the sources (|27|l . The parameter of the 
expansion is g'^[Ui,Vj] [U,V], see Eq. ^^). 



IV. CLASSICAL FIELDS AND LIPATOV 
VERTEX IN THE FIRST ORDER IN [U, V] 

The general formulas for the classical solution in the 
first order in [[/, V] ([Z-/, V]) have the form 

AW<^{x) = ^jd'z({A'l{x)A'%z))A^Xrt{z) 
{AI{x)A'%z))a,j^%\z) 

A'^^^^x) = ijd'z({Al{x)A^^{z))^,j:rt{z) 

{Al{x)A'%z))^^Xr,\z)) (28) 

The Green functions in the background of the Eq. H25|) 
field can be approximated by cluster expansion 

{A,{x)A'^{z))a,^ 

= {A^{x)A''{z))uM + {A^{x)A''{z))vy 

- {A,,{x)A%z))o + 0{[U,VUUy]) (29) 

(and similarly for other sectors) where {A^{x)A'^ {z))q are 
the perturbative propagators (0) and the propagators in 
the background of one shock wave are given in the Ap- 
pendix As a first step we shall discuss the behavior 
of classical fields at the t = —oo boundary. 



Note that while each of the field Ui and Vi satisfies the 
boundary condition Ai(t — > — oo) = pure gauge, their 
sum (|25ll does not. I will demonstrate now that the cor- 
rection (|28|l restores this property so A^"^ +A^^'> is a pure 
gauge (up to [?7, terms). 

We need to prove that D^j_A\}^ — /j, i/ cancels Fik = 
—ig{\Ui,Vk] — i k) a,s t ^ — oo . First, note that the 
contribution from terms [?7i,y*] in T (and [Z^i,V*] in T) 
vanishes at i ^ — oo since these sources are located at 
t > 0. Also, it is easy to see that the Green functions 
{Afi(x)A'^ (y)) interpolating between sectors also vanish 
at this limit (see the explicit expressions in the Appendix 

EI). The only nonzero contribution to A^'' eX t — — oo 
has the form 



Af\x)^gldad(3d:'z^ 



-ictx^ —il3x^ 

' [a — i€){f3 — ie) 



ir\x^\ 



:\0, z^)i[U„Vkh - I k) 



af3s — + ie ' 

k k ,„ 

g{x^\U^U^ + V^V^ ~ ^|[C/„ 14] - * ^ k) (30) 
P± 



PI PI P 

Throughout the paper, we use the notations 

{x^\...\f)= I d^zAxi_\...\z^)f{z^^ 



(x|...|0,/) 



d''z^{x\...\Q,z^)f{z^) (31) 



From Eq. Ip!n|l we obtain 



7(a:^|{/^C/t 
Pi 



V 



Because Fj^ 
p\Fij and therefore 

D,A 



PiP yt _ P^P \abpb 
P± • 

tik in two dimensions pip^ Fjk 



(1) 



(32) 
-i^ j 

j ^ 
(33) 



so the sum Aq+Ai is a pure gauge (up to [U, V]'^ terms). 

It can be demonstrated that if we use background- 
Feynman gauge for calculations in further orders in [U, V] 
parameter we obtain a pure gauge field Wi = Ui+Vi+gEi 
where {idi+g[Ui+Vi, )£'* = 0. Similarly, in the left sector 
the pure gauge field att — — oo will be Wi = Ui + Vi+g£i 
with £ satisfying the equation {idi + g[Ui + Vi, )£J = 0. 

In the lowest order in [U, V] we get (cf. ref. 0,112) 



E. = ix^\U^U^ + V^V^-^\[U..,Vk] 
P\ P\ P\ 



k) 



PI P± 



IT 
Pi 



\[U„Vk]-i^k) (34) 



The explicit expression for the correction Ei in the second 
order in [U, V] is obtained in the Appendix IbI 




(a) 



(b) 



(c) 



(d) 



(e) 




(f) 



FIG. 4: Classical field in the first order. The shaded area represents the linear term 1271 . 



B. Gluon fields in the first order 

To find the amplitude of particle production due to 
the scattering of the two shock waves we should study the 
behavior of gluon fields at t ^ oo. The general expression 
for the gluon fields up to the first order in [U, V] is given 
in Eq. (|28|) and the corresponding diagrams are shown 
in Fig. El 

Let us find the gauge field A^^^ in the x,, x, > sector 
of the space (we use the background- Feynman gauge). 
This field is a sum two terms: due to [C/i,y] and due 
to [t/i,Vfc] — i ^ k. The [?7i,y*] term comes from the 
diagrams shown in Fig. 2t and Since the Green 
functions in the shock-wave background in the forward 
cone a:, , a;, , , > are just the bare propagators |(2J), 



we obtain 



+ ig{xf-t^2iT5{p^)\Q,[U,,V]) (35) 



where the first term in the r.h.s. of this equation comes 
from the diagram in Fig. and the second one from 
Fig. The factor 1/a in the denominators comes from 
the integration over the final point z, from to oo in the 
expression for the [i7i, V"^] part of the linear term l|27|) . 

The term ~ [[7^, Vk] is obtained by integration of the 
Green functions (A1)-(A3) in the diagrams in Fig^jD and 
Fig2^ over the final points z^,,z±_ < 0: 



A,{x) = (36) 

2 J Stt-^ (q — ie){P^ — i€) I aps — p]_ + ap's — + 

k 

+ ix^\ei-a)2nSial3s - pl)d,U——^^ -C/t^^O, [U,,Vk] - i ^ k) 

ap's — p^^ — le 

= 9{x\^-^^^\0, [U,,E^])+g{x\^-^27T6{p^)\0, [U,,P]) 
p'^ + le a 



where we have used the notations (|34() for brevity. The field AI'^ is the sum of Eqs. (|35|l and (|36|l . The field A]} is 
obtained from A,'^ by the replacements and U ^ V. 

Finally, let us calculate the field Ap-* coming from the diagrams in Figs, and^. Using the cluster expansion 
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(I29|l for the Green functions -pjPk and the set of the propagators from Appendix lA II we obtain 



dad(3 



1 



al3s — p-^:+ ie 



dp' 



-U- 



+is 



da' 



-V 



P 



27r(a' — ie) a' 13s — p\ + ie (a — ie)(/3 — ie) a(5s — p\ + i 



2'K{l3'-ie) a(3's-p\ + ie 

ah 

\Q,[U^,Vk\-i^k) 



1 



P 



+i{x^\e{-a)2-K5{a[3s - p\_) 
da' 



dp' 



-IS 



-V 



p- 



27r(/3'-ie) aP's-p\+ie 



-vt- 



27r(a' — ie) a' (3s — -\- ie {a — ie){p — ie) af3s — p]_-\- ' 



ah 



\Q,[U^yk\-i^k) 



2g{x\ 



1 



+ ie 



\Q,E,) + 2^g{x\e{~po)5{p')\0,£^) 



(37) 



The left-sector fileds A^j. are obtained by trivial replacements. 
Let us present the the final set of gauge fields (at x,,a;, > 0) 



(1) 



(1) 



g(x| ^" + ;'^ ' |0, [U^: V% + 2g(:r|^^^^|0, [C/„ E% + ig(xf±^2n8{p^)% [U,,V, + 2£'^]) 



p^ + le 



p'^ + le 



A 



5(^1 ^^^^7^10, n) + 2g{x\ ^\'']' % [F„ E^]) + ig{xfAJl2^5{p^)% [VM + 2E^ 



p^ + «e 
1 



p^ + ze 



A 



2g(x|^— |0,^,) + 2zg(a;|e(-po)2^<5(p^)|0,fO 
+ le 

g{x\^^^^% [U.,V% + 2g(x| ^"~''^ ' |0, - ig(a;|^27rJ(p2)|0, [C/„ V- + 2ii;^]) 

p-^ — ze p'^ — le a 

- g{x\^-^P^\0, H,V]) + 2g{x\^-^^^\Q, [V,,£^]) ~ ^g{xfM2nS{p')\0, [V,„ + 2E'']) 
— le p'' — te p 

2g{x\^^\Q,£i)~2ig{x\e{po)2Ti5{p')\Q,E,) 
— le 



(38) 



where E and £ are given by expressions l|34|l . It is easy The vertex H39|l is transverse: k^R^(k) = kiE'^ + [[/, 
to see that the fields (|38|l satisfy the boundary conditions Vi, E"^] = due to our choice of Ei such that {idi + [Ui 
A{x^ + X, oo) ^ A{x^ +x, oo). The fields Af, for Vi, )E'' = 0. Similarly, fc^7^p(fc) = 



the right sector only (with Ui = Vi = 0) coincide with 
refs. M,M 



It is worth noting that in the case of two weak sources 
Ui ^ {k — k')i and Vi ^ k'^ the vertex reduces to 



C. Particle production and Lipatov vertex 

The particle production is determined by the behavior 
of the fields at t = oo described by the Lipatov vertex 



i?^(fc_L,??) = lim k'^A^ik) 



g[U,,V' + 2E' 



Ps 



(39) 

g[V,U' + 2E'] 



TZf^(k±,T]) = lim k-A^{k) 



2g£^ + ^g[U,,V + 2^] + ^litg[v,^W + 2^ 
as ps 



2 

s 



El^El] +2k' -2fc, ^^'^'^^ 



{k',k-k')^(^ 



(2^' - k)- + (. - - (/? - ^)P., 



Ps 



(k^k') 



i\1 _ Z,'2 



kl 



where the first three terms form a standard Feynman- 
gauge Lipatov vertex and the last term is a longitudinal 
contribution which drops from the particle production 
amplitudes. 
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We get 

R{k]U,V,U,V) (40) 

It is covenient to rewrite this product in terms of the 
transverse part of the Lipatov vertex in the axial P2 ~ 
gauge: 

(41) 

It is easy to see that 7?.pi?^ = £iL' so the Eq. 
reduces to 

R(k]U,V,U,V) (42) 
= /:^(/fci)L^(~fci) e ^Id-^^i-uTVT+uTvn 

The exphcit form of the axial-gauge Lipatov vertex is 

U{UX) = 2gE,~2ig^[Uj+2Ej,V'] (43) 

The apparent asymmertry between U and V in the ex- 
pression (|43|l for Li does not affect the results - alterna- 
tively, one can use the expressions (|43|l with U ^ V,U ^ 
V and the product H42|l will remain the same. 

The longitudial part L can be easily obtained from H39() 
but we will not need it. Note that the dependense on rj 
is governed by the slope of Wilson lines. 



V. kr FACTORIZATION FOR THE DEEP 
INELASTIC SCATTERING FROM THE 
NUCLEUS 



In this section we will reproduce the standard k±- 
factorized result jH EHIE 113 for the number of pro- 
duced particles ^ for the case when the Ja is small (e.g. 
virtual photon) and Jb corresponds to nucleus. 

A weak source Ja produces only one gluon and J^a 
absorbs this gluon so the upper part of the diagram is 
attached to the lower by two-gluon exchange only. For 
the two-gluon exchange, factorization formula H19|) sim- 
plifies to 



DADAJipA)JiPB)e-'^'-^^+'^^^^ (44) 

^ gi/d^ZiVn^i) U"{z±)~iJd'^z^V^''{z^) U"{z±) 

j DADAJ{pa)J{pb) e-*^(-4)+^^(^) 
X d?z^Sz'^ HiV^.V^y - l)(Vl,V, - 1)} 



^ir{{U^Ul,-\){UM-l)} 



1 1 










i 

i 



= < 



• • z 



(a) 



u 



y 



(b) 



FIG. 5: Factorization for the two-gluon exchange. 



where Vz — [ooni + z_\_, — oorii + z±\, cf. Eq. p9l) . This 
formula is easily seen from Fig. Indeed, in the lead- 
ing order in perturbation theory the l.h.s. of Eq. I|44() is 
represented by the diagram in Fig. [S^ where the quark 
line in the external field lAi^ Ui is given by 

ixixMy)) - [dzSiz,)U^{x\hz) ^2UlUz{z\hy)Ul 



(45) 

at x*, < (cf. Eq. 1111) '). On the other hand, the r.h.s. 
of Eq. is represented by the diagram in Fig. [5|d. It 
is easy to see that the factor d'j_ cancels the propagator 
= ^ so the diagram in Fig. ]^ reduces to|3jL. 

If we neglect the evolution, the slope of the Wilson 
lines U can be replaced by pi and the slope of V^s by p2 
(see the discussion in Sec. Ill Bll . The number of particles 
produced in a collision of weak and strong shock waves 
(I24|) is given by the square of the Lipatov vertex 139|l . 
Technically it is convenient to introduce a source Xi{x±) 
for the Lipatov vertex Li ( and a{x±) for C) so 



X e 



For a weak source Vi ^ diV we get gEi ~ —{U{gik + 
2f.)U^Ybyb ^^(^^) ^ 2{x^m,U]p^U^r'\Vt) and 

the above equation can be rewritten as 
R{k;U,V,U,V) = J<fx±(fy±e 



5A^(xi)<5a-(2/i) 
where 



u 



U = Ue ^ 



(47) 



dzl d 



± 



(This expansion is similar to the functional-integral rep- 
resentation of the non-linear evolution as e^'^^Ue~'^'^'^ de- 
veloped in ref. 28]). We need here only the first two 
terms of the expansion in powers of A and a. 

Using the simplified factorization formula (|44|) we get 
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R{k;U,V,U,V) 
1 



X tr{UlU.Ul'l^.'}dl. dl,,tr{{V}V., - 1)(V],V, - 1)} 



- _ 1) y '^^^^i'^^yj-rf^^irf^^l 

- {y\^\z'){euM-U.,UleUyUl)\ 

P± 

X [{x\^\z){eu,ul, ^ U^UieU^ul) 

- ix\4\z'){u.ulr -u.ult^u^ul)]} 

p± 

X e'C^^^-^'-aL ai,,tr{(F;K. - l)(Vl,V, - 1)} (48) 

Without evolution, U = U and V = so the square of 
the Lipatov vertex (|48|l reduces to 



R{kr,U,V) (49) 
' d^z^d^z'^ tvidlV, d\vl +z^ z'} 



m - 1 



X Tr{((fc^|[4,C/]|z^)-((fc^|[4,C/]|z^)) 
P± P± 



X ((z^|[4,C/t]|fc^)-(zl|[4,C/t]|fc^))} 
PI PI 

It is easy to see that in the momentum representation we 
get 

R{k^;U,V) = ^I^Jd^k:td^k^ K^FkLikukr,k) 
X tr{V{k^)V^-k^)}tr{U{k± - k^)U^ {k^ - k^)}{5Q) 
where 

_ kl^ik-k2)l , kl^{k-k,f 



P' 



+ 



(fci - k2)i (51) 



is the gluon-emission part of the BFKL kernel [23 • 

Substituting the result of the integration over central- 
rapidity gluons H49|) into the factorization formula ^ we 
obtain the standard /cr-factorization formula 

ng{kr,v)) = ^I^Jd'kid'k2Kg'^\,^{k,,k2;k^) 

X {A\Ti{V{ki)vH-k2)\A) 

X (B|Tr{t/(fc_L-fci)t/^(fci -fc2)|B) (52) 



where the rapidity dependence comes from the slope of 
the Wilson lines in the r.h.s. This formula was obtained 
in the approximation when we neglect the evolution of 
the shock waves, but it can be proved without this as- 
sumption 2^, ,2a, I26L I27I I . It should be emphasized that 
Eq. (|52|) is valid only in the first order in the [U, V] ex- 
pansion (that is, for the pA scattering). As we shall see 
below, it is not valid beyond the first order. 



VI. EFFECTIVE ACTION 

In this section we get the first-order effective action 
for the double functional integral and check that the 
Lipatov vertex H43(l serves as a "splitting function" for 
the non- linear evolution equation |30ll3ll | 



{x-z)\{y-z)\ 



x(Tr{C/,C/]}Tr{t/,C/n - N,!Yy{UJJI}) 



(53) 



The effective action is defined by the functional integral 
(Pljl without A^i{p)A^{~■p) insertion: 

e'S.«'^^-y'U.y) ^ JdADA e^p\^~iS{A) (54) 

+ I jd^z{yT[oo,T.„,oo\l + (oo,.F„,cx))JZif ) 

+ iS{A) - ild^z{v^'[oo,F,,,oo]1 + (oo,F„,oo)^C/r)} 



Performing the shift A ^ A + Af\ A^ A + Af^ we 
get (cf. Eq. 123 

is^ii{u,v.uy) - 



DADAey^^^i j d^zi^i^A^V^^A" 



1 



-A'^D^^A'' -T^A^ + Tf^A^)} 



(55) 



The effective action in the first nontrivial order is given 
by the integration of linear terms with the appropriate 
Green functions 

iS,fiiU,V,U,V) = 

\ j d'zd'z'[ - T;{z){A^'^{z)A'^\z'))Tt{z') 
+ 2T;{z){A^'^{z)A''\z'))%\z') 

-t^{z){a^'^{z)A''\z'))t!;{z')] 

The first term here can be taken from [13. |2^ : 



(56) 



a,A?7 j d^zi^E^E"' + [C/„ V + 2E']^%, W + 2E']j 

where Arj — rji — r/2 is the difference of rapisdities of 
the slopes of U and V Wilson lines. The third term 
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^is^ii{u,v.uy) 



X / DAcxp 



S^g{U,V) in Eq. is obtained from Eq. by the field A coupled to Lipatov vertex 
usual replacements U ^U, F V 

Let us calculate the second term beginning with the 
contribution ~ [J7i, V']. The Green function in 

the x^,x, > sector is simply the perturbative propaga- 
tor (O so one obtains 



(62) 



'9^ I 7S^(0, V.r\2^S{a(3s - pl)d{a)\0, [U\ V^^) 



1 



' da 



(57) 



The integral H57() is formally divegrent at both small and 
large a. This divergence occurs because we've put the 
slopes of Wilson lines ei and 62 (see Eq. (|2U) to be pi 
and p2. If we keep the slopes off the light cone, we get 
(seeref. H): 



dadp 



27r2 (a + e-''i/3)(/3 + e''2a) 
X (0, [U\V,r\27T6{al3s ~ pl)e{a)\0, [U\V^'') 

= - 2a,A7j{[U\V,r\^[U\V^'') 



(58) 



The contributions of [Ui, Vk] ~ i ^ k terms can be calcu- 
lated in a similar manner by integration of these terms 
with appropriate Green functions in Eq. ]b<i\i . After some 
algebra, one obtains 

-2a,A77{ J d^z^rE,{z^) + {W - V^,£T\^\[U\V,T) 

+ {[W,V,r\^W-V,,Er)] (59) 

The sum of Eqs. (|58|l and H59|l give the second term in 
the effective action 

-a,AT] jd^z(2£'''Ef + [U,,V' + 2£']^%, W + 2E'] 
+ [V^M' + £']^[U^,V' + E']^ (60) 

It is easy to see that the total effective ac- 
tion s,s{u,VMy) = s':^{u,v) + S''^{U,VM:V) + 

S'^g (Z//, V) can be represented as 

iS'cff(C/,y,W,V) = - iUX + iU^V- (61) 
+^A7/ j d^zl^L^L""' - 2£"L^ + £'^'4") 

where and Ci are given by Eq. (|43|l . Note that 
the effective action is a square of Lipatov vertex (|39|l : 

For the future applications we will rewrite the effective 
action H61|) as a Gaussian integration over the auxiliary 



A single field Ai for both Li and Ci reflects the fact that 
gauge fields and coincide at t = 00. 

Let us prove now that the effective action 1)62(1 agrees 
with the non-linear evolution equation. To find the evolu- 
tion of the dipole U^U^, we need to consider the effective 
action for the weak source V . As we mentioned above, at 
smah Vi - d^V one has gEi ~ -{U{gik + ^)W)''''V^ 

so L,{x±) ~ 2{xj^\[p-,U]p''U'')''''\Vj^) and Eq. ^ can 
be rewritten as 

J DAD A exp I - iS{A) + iS{A) (63) 
- I J d^z^Vf [c5o,J^„,c3o]^ d^z^{oo,J^„,oo)'^U^' 
+ ^ fd^z^V,°''[oo,F„,oo]l + i fd^z^{oo,F„,oo)lU^'^ 



DA exp <^ d^z 



( 



AfA" 
asArj 



)} 



where 



2-2fA -2-^{U^\U) 

U ^ Ue ^ 

U ^ Ue ^ 



(64) 



Again, we need here only the terms up to the second order 
inA. 

To find the evolution of the dipole UxU^ we should 
expand Eq. (|53|l in powers of Vi{Vi) and use the formula 

(U^Ul)"^^ = Ppi3/J'<iz.[oo,F..,oo],^ YVe get 

DXe-^-f""''^^'^" U^UlUyUl (65) 

Performing the Gaussian integration over A one obtains 
after some algebra 



(66) 



tr{U.UlUyUlr + ^fd^z, 



47r2 J - '^-(x-z)i(z-y)i 
X {tr{UlUxUlU,}tr{UlU,UlUy} - N,tv{UlUxUlUy}Y^ 

which is the non-linear evolution equation for the double 
functional integral for the cross section |2^ [s^ • 



VII. CLASSICAL FIELDS AND LIPATOV 
VERTEX IN THE SECOND ORDER 

If we neglect the evolution, the classical sources Ui and 
hii coincide (and similarly 1^ = Vi) so the corresponding 
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^ >^ ^ 



(a) 



(b) 



(c) 



FIG. 6; Retarded classical field in the first order 



fields in the right and left sectors coincide and are deter- 
mined by the integrals of the retarded Green functions 
with the T = T sources. 



A. First-order gluon field in the x'^^A^ — gauge 

In the first order in [U, V] , the classical fields (in 
the background gauge) reduce to 



+ V,{x^)0{x,)0{~x,) + 2g{x\ 



1 



|0,i?.) 



+ lepo 



{a — ie) ^ 
p2 + iepQ 



7(^1 



p2 + iepQ 



io, [c/.,n) 



2g{x\ 



1 



(67) 



We see that the fields outside the forward cone are piece- 
wise pure gauge: 

A^'^ = w^{x^)e{-x,)e{-x.)+ u^{xi_)e{~x,)e{x,) 
+ v,{x^)e(x,)e(-x.) 



(1) 



1 



ig5{x,)6{-x,){xi_\^\[U,,E']), 



A 



(1) 



ig6{x,)6{~x,){xi_\ — \[V,,E% (68) 
PI 



while the field in the forward sector a;* , a:, > is deter- 
mined by the Lipatov vertex H39|l . 



p^ + lepo 



(69) 



Summarising, the classical field H67|l can be represented 



as 



where 



4'^ = A^, + {x\ \ !<)) (70) 



A, = e{-x,)e{x,)u^ixi_) (71) 

+ 9{x,)0i-x,)V,ix^) + e{~x,)0{~x,)W,{x^) 



is a trivial part corresponding to a piece- wise pure-gauge 
field and 



^ ^ s \a + ie a — It 



s \p + le p + le J 



s \(3 + ie 



f3 + ie 



(72) 



describes the non-trivial part related to the gluon emis- 
sion. 

The axial-gauge Lipatov vertex is given by first line in 
Eq. (ESI: 

L,{U,V) = 2gE,-2ig^jUj+2E,,V'] (73) 

Following ref. it is instructive to represent the fields 
lIHTjl in the gauge xjl^l^ = x^,A, -\- x,A^, = 0. In this 
gauge one obtains at , x, > 



Al{xi\,x±_) = -J udu xPFpf,{uxi\,x±), 

nOO 

A-^{x\i,x±) = du x'^^Fpf,{ux\\,x±) (74) 



where = -I- ^x,p2 (so that x = x^^ + x±). The 

limit of integration oo in the above expressions was cho- 
sen to satisfy our boundary condition (no pure gauge 
fields at t = oo 

The fields ouside the forward cone (|68|) trivially satisfy 

the xjiA^^ — gauge condition. The fields in the forward 
cone x, , x, > are obtained by integrating f 's corre- 
sponding to the fileds in the bF gauge (|67|l . From the 
Eq. ((73|l we get: 

A,{x) = [d^z^e{-x\ + {x-z)\) 

X {2[U,,V'] + [U,-V,,E']) (75) 

A^V(^) ^ _ j rf2^^^(_ 2 ^ _ ^^-^2^^ 

X {2[U,,V'] + [U,-V^,E% 
A,{x) = -^fd^z^5{x\~{x~z)l)E, 



j2_ Of 



+ 9 / d'z^B{-x\ + (x - zY^)[x^\^\zx_)\V, ^ Vj.E'] 

Similarly to ref. 0, the fields x^A'-J-\ x,Ai^^ and ^ 
are boost invariant. However, as we mentioned in the 
footnote, the fields (|75|l differ from those in ref. ^] due 
to a different boundary condition. 



B. Gluon field and Lipatov vertex in the second 
order in [U, V] 

In the next order the classical field ^^^^ is given by 
diagrams in Fig. [7| calculated in the Appendix C. The 



(b) ^ ^ (c) ^ ^ 



(d) ^ ^ 



.7^ 

(e) 



(f) 



FIG. 7: Classical field in the second order in [[/, V\ 



result of the calculation is given by the sum of the piece- 
wise pure gauge field and the field of the gluon emission 
described by the second-order Lipatov vertex represented 
by two terms coming from the diagrams in Fig. El a-nd 
Fig.d 



^ 0{~x^)9{x,)Ui{xj_) + e{x^)e{^x,)V^{xA 



(76) 



+ {x\^^. — 
p'' + lepo 



The first part of the Lipatov vertex coming from the di- 
agrams in Fig. Eland Fig. [TJi has the form 



(1)/ 



(77) 



The second-order term coming from the diagrams in Fig. 
[71;-f is given by 

^ \as ps/fi 

(k',k- k')A^ + !^ -iy/g . , 
- ' ' ^-R\{k'^)Rl{k^-k'^) 



k'\{k-k')\ 



k'l V f3s 



V PS as 



X gE,\'{k^ - fcl) 



(79) 



/3s as / 
where we use the notations 



) {k^\[U,,vr'^\[U,+V„E^'') 
/ ^^ p^ 



/fell 



K={k- 2k')j_ + —-{k, k - 2k') 
aps 



(80) 



where the notations are 



{R['lr 



2gE^^ + A^g{i^^U) — U^r''[V,,E^]' 
P± 

+ Azg{{d^V)^V^r'[U.,E'f 
PI 

1 



= - g[U,, Vr - 2g((dlV) — V^r'[U„ E 



111 b 



2gpl{U^U^r'[V,,E^]'' 



g[U,,Vr - 2g{idlU)^U^r''[V,,Ef 

P± 

2gpl{V^V^r''[U.,Ef (78) 



g = k'l{k^k')l-[^^^^ + {k',k-k')^)y (81) 



(82) 



i?^(fcl) = 2£;^e(^l) 

ik,k')± + ^-iVg f2pi p . ^ 2p2 , 
+ lil + ^^^('-^^.c 

and i?[ ^ = r\^ + r\}_, r\ = -^2+ -^2—' 

With [U,Vf accuracy R^^^ and i?^^^ in Eq. ^ can be 

simplified to 

{R^Pr = ^ g[U,,Vr + '^9Pl{u\u^T\V,,E'f 

P± 

{Ri'^r = g[U,,Vr + 2gpliV^V^r'[U^,Ef {83) 

P± 
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The second-order Lipatov vertex is the the sum of Eqs. 
(|77jl and (O at the mass shell P = 0. At the first sight, 
it looks like the expression H79(l is divergent a,t k'j_ \\ k± 
since Q = k^k'^{l — cos^ 9). This coUinear divergence is 
however purely longitudinal and therefore can be elimi- 
nated by proper gauge transformation. To see that, let 
us write the Lipatov vertex in the axial light-like gauge 
pi^A^ = (|41|) . As we mentioned above, only the first 
transverse term in r.h.s. of Eq. (|41|) is essential since 
P2 term does not contribute to the square of the Lipatov 
vertex. For this transverse part we obtain 



(84) 



-(1) 



(k^) = R^\k^) + ^{RYi{k^)+RY>{k^)) 



2k j 



if) ^(fc^) 



abc 



(PkY 



47r^ 



fcl k' 



{k,k 



l\2 



(85) 




^) {2k' ~ kyL\'^\kY)L^^^'\k^ ~ kY) 



i(fc,-i^fc;)L«''^(fc^)Lr)^(fc^-fcl) 



abc 



jfk^ 



2[U,+V,,E^\\k'^)Et{k^-k'^) 



2^[U, - 2Eiy]\kY)E'^Ak^ - fcl) 



,(fc 



-|[c/, + y„i?Y) 



Note that since the square bracket in the r.h.s. of the 
above equation vanishes at k'^ \\ k±, the coUinear diver- 
grnce is absent. The second-order Lipatov vertex H84|l is 
the main technical result of this paper. 



VIII. CONCLUSIONS AND OUTLOOK 

Let us summarize the progress towards the solution of 
the main problem - the particles/fields produced in the 
collision of two shock waves. The Yang-Mills equations 
with sources U and V describe the two shock waves corre- 
sponding to the colliding hadrons. The expansion of the 
classical fields in commutators [U, V] has the advantage 
of being "symmetric" in contrast to the usual expansion 
in powers of the strength of one of the sources. We have 
calculated the second nontrivial term of the expansion. 
This term is relevant for the description of dA scatter- 
ing, similar to the first term ~ [U, V] describing the pA 
collisions. 



(a) 





(b) 





(c) 



FIG. 8: Classical field for W / i7, V / V, 

Note that while the first-order field given by Eq. H7()|l 
(or Eq. (EHJ) is real, the second-order field has an imag- 
inary part given by the second term in braces in Eq. 
(|79|l . The real part of the second-order term is given 
by Eq. H77|l plus the first expression in braces in Eq. 
H79|l represented by the product of first-order Lipatov 
vertices. I think that this univeral structure will survive 
to the higher orders of the commutator expansion. Un- 
fortunately, the explicit form of the imaginary part of the 
field (second term in the r.h.s of the Eq. ijSSl) does not 
suggest any idea how this expression may look in higher 
orders in [U, V] expansion. Technically, the relative sim- 
plicity of the real part is a consequence of its relation to 
the leading log approximation (LLA). If we consider the 
general case U ^ U , V ^ V , the second-order classical 
field would contain A77 just like the effective action (|61|l . 
Indeed, if we calculate the field in the right sector, 
the typical expression 



du 



{u + ie)Z{u) 



[U,V]{k')[U,V]{k-k') 



(see Eq. HC16|) ') would be replaced by 



du 



{u + ie)Z{u) 



du 



3 (u + ie)Z{u) 
du 

[u -t- ie)Z{u) 



[U,V]{k')[U,V]{k - k') 

'U,V]{k')[U,V](k-k') + 
[U,V]{k')[U,V]{k~k') (86) 



where the three terms correspond to the contributions 
of diagram in Fig. |SjL,b, and c, respectively. The inte- 
gration over u is regularized by the width of the shock 
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wave (cf-li^lsll)) and only the real part is survives - the 
imaginary part exceeds the accuracy of the LLA. 

If we consider the amplitude rather than the cross sec- 
tion, we take only one set of fields (to the right of the cut) 
and impose the usual Feynman boundary conditions. In 
this case the classical field is the sum of logarithms of 
the type (a^ lns)"[J7, F]" (cf. 0). The imaginary part 
calculated above may be related to an old idea due to 
Lipatov that one can unitarize the BFKL pomeron if one 
finds the proper iir or —in to each In s in the LLA ap- 
proximation. Indeed, both these imaginary parts come 
from one source - causality: the itt's in the amplitude 
come from the dispersion relations based on causality, 
while itt's in the classical field H79|l come from retarded 
propagators. 
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while at , < the propagators are 

{A,ix)AAy)) = UAx\-^\y)Ul (A2) 

p'^ — le 

{A^{x)A,{y)) = -UM^^O^,{WU)^^\v)Ul 
where 

Op.(C/) = JdzSi^z,) \z){2ag^,U (A3) 

s as"^ ) 

2. Retarded propagators 

First, let us present the retarded propagator in 
the background of the shock wave 0(a;*)rii(a;j_) + 
0{—x^,)Ai{x±) where fti = ^li^ift'' and A^ = AidiA^ are 
the pure gauge fields. This propagator can be obtained 
from Eqs. HA1|I . HA2|I by setting U ^ U ^ il^A, taking 
appropriate combinations, and rotating by the matrix f2: 

{A{x)A{y)),,, = 6{x,)e{y,)^M \ \y)nl 

^ Q{-x:)d{~y.) A,M \y)^l (A4) 

+ g(a;,)g(~yO»x(x| \ Op.(»U) \ |y)At 
p^ + lepo p'^ + tepo 



APPENDIX A: GREEN FUNCTIONS IN A 
SHOCK- WAVE BACKGROUND 



1. Feynman rules for cross sections in a shock- wave 
background 



Let us present the set of the bF-gauge propagators in 
the background of a shock-wave field Ui9{—x,,), Ui9{—x^.) 

At x,, > all propagators are bare, see egn. 
At > 0, < we get 



{A,{x)AAy)) = {x\^]--0,AU)^^\y)UliAl) 

(^^(x)A(2/)) - {x\^^0^,{U)^-\y)Ul 
{Af,ix)A^{y)) 



{A^{x)A,{y)) 



^{x\27^Sip^)0ipo)O^,iU)^^\y)Ul 
p^ + te " 



1 



iix\27rSip^)9{^Po)0^,iU)- -\y)Ul 

p^ — le " 



a. Cluster expansion 

The background field in our calculations is the trial 
configuration Ai = d{—x^)Ui + 9[—x,)Vi. Since we are 
expanding in powers of commutators \U, V] , the adequate 
procedure for the propagator in the Ai background is the 
cluster expansion (|29|l : 



(A5) 



p2 + iepoJ u+v 
/ 1 \ / 1 \ 1 ^ 

\P^ + iepoJu \P'^ + iepQ/v p'^ + iepo 

where dots stand for the second and higher terms of clus- 
ter expansion. Most often, the first term (|A5I) is suffi- 
cient. In several cases when we need the second term, 
the fillowing trick helps. 

Let us add and subtract Ei to our trial configuration 
so it takes the form Ai = A^ — Ei9{—Xt,)9{—x,) where Ai 
is the piece- wise pure-gauge field given by Eq. (|7TTl . see 
Fig. IHI With [J7, accuracy, the propagator in the Ai 
background takes the form: 



1 



1 



1 



P2 + iepo p2 



lepo 



p2 



lepo 



-{p\E.,} 



1 



P2 + iepo 



(A6) 
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X_ 



FIG. 9: Piece-wise pure-gauge field Ai 



where we can replace (P^ + iepo)^^ in the second term 
in r.h.s. by the first term in cluster expansion (|A5I) . The 
remaining first term in the r.h.s. of eq. (jA6|l is calculated 
below. 



b. Propagator in the piece-wise pure-gauge field. 

The retarded propagator {Af^{x)Ai^{y))yct in the back- 
ground of the piece-wise pure gauge configuration shown 
in Fig. I^can be obtained by "squaring" of the propa- 
gator in the background of one shock wave HA4(I . In the 
region a;* > 0, 2;, < and y*,?;, < it is given by Eq, 
^4|) with appropriate substitutions: 

{A^{x)A,{y)),,t (A7) 

+ lepa -I- lepo " 



where 



(A8) 



dzS{-z,) \z){2ag^,V^W + ^{p^.d^V^W) 
as'' J 



Here the last term ^ [Ui , V''] (additional in comparison 
to Eq. IjASp ) is due to the source contribution to the 
second variational derivative of the action 



5^S 



(A9) 



SA-{x)SA''^{y) 

+ i[V,)W + i<5(-x.)(a, + t[Ui)V'}] ''"six - y) 
(is J J 

The propagator in the region a;, < 0, a;, > and 
2/*,?/, < is similar to Eq. HA7|I 



{AAx)AAy)),^t 



(AlO) 



1 



+ iepQ 



1 



p^ + iepo 



\y)wi 



where 

O^AU^W) (All) 
''dz5{-^z,) \z){2Pg^,U^W +jip2,dAU^W) 



{dl{u^w)-^i[v,,W]w)}{z\ 



Finally, the propagator at , a;, > and y* , y. < 
has the form: 

{A^ix)A,4y)),,t (A12) 

= t{x\^^. — \o^AU)^^. — 0«,([/V) 
^ 'p^+tepQ I '^"^ 'p^ +iepa ' 

p^ + zepo > V + «epo 



APPENDIX B: PURE GAUGE FIELD E IN THE 
SECOND ORDER 

From F^,^,(l}^ + + S^) = we get 

(a, - i{W,, -t^j = z([C/„ Vj] -1^3- [E^.E^]) 

(Bl) 

If we choose the {di — i[Wi, )E^ = condition the above 
equation reduces to the recursion formula 

E^ = -g{x^\W^W^\^\U,,Vk]'^i^k~g[E,,Ekf) 
PI 

m 

It is convenient to introduce complex coordinates in the 
2-dimensional plane: z = zi -\- iz2, z = zi ~ iz2 and Q — 
Qi -\- iQ2, Q = Qi — iQ2 for arbitrary vector Q. In these 
notations the recursion formula (|B2|I simplifies to 

E" ^ ~g{xAW^W^^''Ki2~Ei2) (B3) 
P 

where K12 = [Ui,V2] - [U2,Vi] and E12 = [Ei,E2]. In 
the leading order W^W^ can be approximated by cluster 

expansion: W^W^ ~ ("^)(i) "^^^^^ (^) ^ ^^^^ 
V^V"^ — 4 and therefore we get Eq. (|34|l . In the second 
order we need one more term of the cluster expansion: 

W-W^ ~ ( J-) + {U-U^ - -)p{V-V^ - -) 
P ^P{i)^ P P P P 

+ iu\u^ - l)p{V^V^ - i) - (^)e^'^ (^) 

P P P P ^Mi)^ 



so the second-order expression for E is 
1 



ig[x^ 



PJH 



"'[ElE^]'^) 



(B4) 



ig(x±\iU-U^ ~ -)p{V-V^ ^-) + U ^ V^r^ifJa) 
P P P P 
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FIG. 10: Classical field in the backwards cone 
Similarly, for the E component we get 



(B5) 



P P P P 
The corresponding formula for the pure gauge field £ 
in the left sector is obtained by the trivial replacements 
U andV ^V. 



APPENDIX C: CLASSICAL FIELDS IN THE 
SECOND ORDER 

1. The fields at x,,x, < 

Since all the Green functions in our expansion are re- 
tarded, the only second-order contribution order the clas- 
sical field A^^) is comes from the D^Fik part of the linear 
term shown in Fig. ^| At < the gluons in Fig. 

[m propagate in the external field Ui ^Vi. It is conve- 
nient to add (and subtract later) the external field Ei. 
The contribution of the diagram in Fig. 1101 gives then 



i I d}z{x\ 



{P - E)^g,k + 2iG, 



-\zr'D^L%{z) 



k k 

= -{x^\W^W^\'^'-L\^) + 2{x^\W^W^\'^''[E,,Euf) 
PI PI 

+ (x^lW^W^r'^iUk.E"]'^) 
P± 

Since each of the two legs in the diagram in Fig. [7^ 
represents Ei with our accuracy, the Fig. contribution 
can be reduced to 

{x^lW^W^l'^hiD.Ek - I ^ k^E^^f) 
P± 

k 

- ix^lW^W^-'lE^^Ek]') 
PI 

Combining these two terms we get (at , x, , 0) 



Pa 

k 

+ (x^lW^W^l'^'iE^Ek] 
PI 



El{xi_) (CI) 



up to the terms '-^ \U,VY', see Eq. HB4|I . Also, it is 
easy to see that the longitudinal components A, and A.^, 
vanish at x^,x, < 0. 




(a) (b) 

FIG. 11: Classical field iat a:* > 0, a;. < 

2. The fields at x, > 0,a;. < 

First, we note that there are two types of diagrams 
shown in Fig. 1111 with the three-gluon vertex in the 
z» > 0, z, < quadrant and in the 2;*, < quadrant. 
The contribution of the first type (see Fig. Illb ) vanishes 
because the only non-zero component of the first-order 



field J^^' in this case is A^^^ such that 13* 



(1) _ 



0, see 



[(1) 
the Eq. 

Next we calculate the diagram in Fig. Eb- As in the 
previous case, the gluon legs are attached only to the 
D'^Fik part of the linear term H27|l . The gluons in Fig. 
[TTb propagate in the external field Ui9{—x^.) + Vi9{—x,). 
The propagator in this background is given by the cluster 
expansion IjASp or, if one needs the [U, V]'^ accuracy, by 
Eq. ((X6|) and formulas l|J7|l - lPl2|l . 

Let us start with the transverse component of the field 
Ai . If the three-gluon vertex is integrated over only the 

, z, < quadrant one can demonstrate that similarly 
to the x^,,x, < case, the contribution of the diagram 
in Fig. [TTb reduces to 



d^z0{-z,)e{-z,) 



(C2) 



1 ^, 



{x\^P'\zr''[Ll,{z^)-[E,,E,]\z^) 



where Pi = + gA^, Ai = UiOi-x^) + V,e{-x,) + 
Ei9{~x^,)0(—x,). Using the Green function in the the 
two- shock- wave background (|A10|I . we see that the r.h.s. 



1(2) 



0. 



of Eq. Id vanish so A\ ' = A 

It is easy to see that = at a;* > 0, x, < so we 
are left with A, only. Again, since the only contribution 
from the three-gluon vertex comes from the z* , z, < 
cone, it can be demonstrated that 



Ai{x) = ~2ijd^ze{-z^)e{-z,) 



(C3) 



X {x\j^F.~P^\zr>'{Uk{z±) - [E.,E,]{z^) 



Substituting the explicit form of the Green function 



18 



b 



(|AllA2p one obtains 

-2i J d^ze{-z,)e{-z,){x\V^V^F„W^W^zr'' 

= 2{x\V^V^^—r''0,[U,,Ef)6i-x,.) 

a — te s 

Thus, the only non- vanishing component of the classical 
filed in the x» > 0, < region is 

A:ix) -2i5i-x,){x^\V^V^''\U,,E'f) (C4) 
s 

Similarly, at a;* < 0, x, > 

A:{x) = -2iSi-x,)ix^\U^U^-'^[V,,E']'') (C5) 
s p^j_ 



and the last term 



s{a + ie) 



([{/„ V']U^U^Vj)''''E^'' (C8) 



is actually a part of the second-order contribution coming 
from the [C/i,y*] term in the Green function (see Eq. 



The remaining part of the Lipatov vertex (coming from 
the diagram Fig. with z outside the backward cone) 
will have the same structure H77|) with different R± , Ri , 
and i?2- Note that our first part of Lipatov vertex H77() 
satisfies the condition 

P^^rW (p) = + r['^ + r['} + R^l + i?W = (C9) 

(recall that {idi + g[Ui + Vi,)E'' = 0). 



3. The fields in tlie forward cone a;,, a;, > 

a. The three-gluon vertex in the backward cone 

At first, we consider the contribution from Fig. [7^ 
where the three-gluon vertex integrated over z* , z, < 
quadrant. This contribution is similar to the one we con- 
sidered in the previous section so we can start with the 
expressions (jC2|l and (|C3p . Using the Green functions 
in the A background given by Eq. (|A6() and formulas 
(|X7|) - (|A12p . one obtains 



(C6) 



J p-^ + lepo L 

^(Pit^^ + ,,dM - Elt^] ^[/t|-&0, [K„ Ef) 
s\p — ie a + ie J p\^ 



s\j3 — ie 



b. First part of the Lipatov vertex 

The sum of all the contributions calculated up to now 
(which includes everything but the terms with three- 
gluon vertex outside the backward cone) can be rewitten 
in the form of Eq. l(75|l 



where r'}}^ {p) is given by Eq. lf77|) 



^ s \ p + le p — le 



s \ a + ie a ^ ie 



(C7) 



c. The three-gluon vertex outside the backward cone 

Here we must calculate the diagram in Fig. with 
three-gluon vertex outside the backward cone a;*, a;, < 0. 
With our accuracy, each of the two legs in Fig. |2b can 
be represented by the field 



A, 



(^1 — 



lepo 



(CIO) 



so we get 



= i J d^z{x\{P^g^„ + 2iF^^ + iepo)-'\zr'' 
X [A('^^,D^A^^'^ ~2DpAi'^]\z) (Cll) 

where we have used the gauge condition D^^A^}'' — . 

It is easy to see that the term ~ G^a in the above 
equation can be dropped. Indeed, since the point z lies 
outside the backward cone, the only non-vanishing con- 
tribution proportional to, say, Fti can come from the 
quadrant < 0, z, > where the only surviving compo- 
nent of the field is = -i6{x,){x±\^\[V^, E']). (see 

eq. (|68|l '). In addition, D,Ai^'' = and therefore all 
possible terms ~ G^^ in I'-h.s. of Eq. (|C11|I vanish and 
therefore 



Al^' = e{-x,)ei~x,)w,ixi_)+ e{~x,)e{x,)u,{x^) a(;}\x) = 
+ e{x.)e{-x,mx^) + {x\ \ + <5i?f) 

+ lepn ^ ^ 



id z{x\ 



P2 -I- iepQ 



(C12) 

\zr''[A^^>,D^Ai^^ ~2DM'^]\z) 



For the same reason, the Green function „-> }■ — in r.h.s. 
of Eq. (jC12|l can be replaced by bare propagator ^^^^.^^^ . 
Indeed, these expressions differ only outside the forward 
cone which means either z* < 0, z, > or z* > 0, z, < 
quadrants (recall that we exclude the backward cone 
z*,z, < 0). Consider the contribution to r.h.s. ljC12|) 
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coining from < 0, z, > quadrant. The only nonzero 
component of the field in this quadrant is A, (see 
above) and since D^A, = the r.h.s. of Eq. ljC12|) 
vanishes. 



We get 



id z{x\ 



1 



+ iepo 



(C13) 

lzr''[i(l)^9^A«-2a.AW]''(z) 



The Lipatov vertex is represented by the two terms in 
square brackets. We will calculate them in turn. 

The contribution to Lipatov vertex from the first term 



4a)(^) = ^ld'ze^'nA^''>-,d^A[}^{z) (C14) 



First, let us calculate the part of this integral coming 
from the product of two r'"^^ terms. We have 



167r4 (fc'2 + iek'Q)[{k - k'Y + ie{k - fc')o] 
I i'd^k', (l-2w)(api+/3p2)M + (fc-2A:')^ 



2/ 167r3 



du 



(fc' — ku)"^ — {k^ + ieko)uu 



The remaining third part of the second-order term is 
R%{k) = -2*/d4^e''=^[Q(i>,9.QW](z) 



be 



•2ij r°° 



d^k 



du 



- 4. 



{k - k')l{u + ie) 



X R[l\k'^)4l\k^-k'^) 



+ 2{hg^, - t ^ j)Rll\k'^)R'^J{k^ - k'^) 



(C17) 



Combining Eqs. ifUTsI) . ((UTe)! . and ifHTT)! we get the 

total second-order term 



4^) \k) 



, d^fc^ r du 



((1 -2u)(api +/?p2)^ 



16^3 7_^Z(m) 
+ (fc-2fc')^)i?W^fc^)i?«^'(fc^-fcl) 



s(/3u -I- ie) s{au + ie) 
u{(^Pl-/3p2-k^ + 2k'^)^ \ n) (1) 



+ 2(fc,9„ - z <-> j)R^\k'^)RfJ{k^ - fcl) 



s V /3 + ie a + ie 



(C18) 



( f ^ + — Vfc, /c - 2fc')± + (fc - 2fc')^| (C15) The integral over u yields Eq. 
lyps as/ } u 



Note that k 



0. 

?(2) 



The second part of -R^(i) is easlily calculated using 
formulas from Appendix D with the result 



-^[i(12)(^) 



abc 



d'^k' du 



Wtt"^ J Z{u) \s{f3u + ie) s{au + ie) 
a(api-/3p2-fc±+2fc^)^ x (1) (1) 



{k - k')]_{u + ie) 

4P2M 
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PI 
1 



s{a + ie) 



m,V']V—V^)^''[Uj,E^]'' 



(C16) 



where Z{u) — [k' — ku)\ — [k"^ + ieka)uu and R\ 



(1) 



APPENDIX D: GLUON EMISSION BY TWO 
WILSON LINES IN THE SHOCK- WAVE 
BACKGROUND 

1. Classical field induced by a single Wilson line in 
the shock-wave background 

In the applications it is sometimes convenient to have 
the result for the classical field and the Lipatov vertex 
in a "non-symmetric" form explicitly expanded over the 
strength of the weak source ~ Ti{p{x±)U{x±y\. This 
expansion corresponds to the diagrams with a gluon pro- 
duction by Wilson lines || pi in the background of the 
shock wave. In this case it more convenient to present 
the results for the covariant gauge shock field A, ~ 5{x^,). 
(The rotation to the pure-gauge field Ui9{—Xt,) is trivial). 

In the first order the classical field is given by the two 
diagrams in Fig^]As in Sect. IVIII we consider here the 
case of the causal classical field corresponding to U = U 
which is the case when we neglect the evolution. Note it 
is not difficult to restore the result tor U U similarly 
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FIG. 12: Emission of a gluon by a Wilson line in the shock- 
wave background. 



to Eq. H38() - roughly speaking, one should replace {p^ + 
iepo)-^U by {p^ + ie)-^U + 2m6{p^)6{po)U. 

The expression for the classical field produced by one 
Wilson-line source can be read from the (retarded) prop- 
agator in a shock- wave background (|A4(I . At > one 
gets 



i{A^^{k)[(X)pl,-copl]z)^at 
R 



1 



fc2 -t- iefco 
^ PI ^ Pi_ 



(Dl) 



-U 



PI J 



s{P + ie) I as 

The emission of gluon by the c.c. Wilson line V^t — 
[— cx)pi, oopi]^ differs from Eq. (|D1|I by sign and the re- 
placement t''U ^UH'': 



i(A^(fc)[-oopi,oopi]z)i.ct 
R 



Kik), (D2) 



{/)_Mai[/^||0,z^r¥C/t 
/ as v : ■) 



s{(3 -\- ie) I as 

The (transverse) Lipatov vertices in the light-like 
gauge are obtained from Eq. (|lT|l : 



4'i(fcx;z±) = 2(fci|[4,C/]|0,Zi)"''C/,<^ (D3) 



= -2(fc^|[4,C/]|0,z^)''Wt 
P^ 

Note that the fields in this Section are presented in the 
bF gauge in the background of one shock wave U which 
differs from the bF gauge for the background field Ui -f- Vi 
used in the bulk of the paper. However, the final result 



■Pi 



(|D3|I for the Lipatov vertex Li corresponds to the P2 — 
gauge and therefore agrees with Eq. (|73|l . Indeed, in 
Sect. IVII it was shown that at small Vi Eq. (|73() reduces 
to 

LW(fc) = 2(fc|[4,C/]/C/tr''0,^,'') (D4) 



which agrees with Eq. ljD3|) if one uses the for- 
mula Vx^y = Pexpig dx^Vi for V{z^) ^ [oopi + 
z±, — oopi + z±] as in Sect. IVII 



2. Classical field and the Lipatov vertex due to the 
two Wilson lines 



In the second order, the field due to two Wilson lines 
is given by the diagrams shown in Fig. 1131 These dia- 
grams are calculated using the retarded Green function 
(|A4p integrated with the three-gluon vertex. The calcu- 
lation is similar to that of Appendix C and the result has 
the form (the details of the calculation will be published 
elsewhere) : 



{A'^{k)[(X), -00]^[-00, OOpi]^')ret 



fc2 -|- iekg 



87r2 



^{k^5^, - k-'Sl ~ ^5«")i?^^(fcl)i?-^,^(fc^ - k'J 



2i 



(El _ El] 

\as Ps/fi 



{k',k-k')^ + !^-i^ 



k'\{k-k')\ 

X [R\y(k'^\ z±_)R^yt {k±_ - k'j_; z^ - i?2y (^1; zx_) 
X R\y,(k^-k'^wjS)-R\v{k';z^) 
X r^t - fcl; ^1) + r\{k'- z^)R^yf (fc - fc'; z') 

- Hg^i - 2 



'''^)[rvHk';z±)VvAk±-k'r,z'^) 



rv{k';z±))r'^^{k± - k'^;z'^) 



where is given by Eq. l|5n|l . and 



(D5) 



= {k'A[pi,u]^\z^r'u.t^ 

p± 

R^2v{k'±;z^) = -^(k'^\dlu^\z^r'u,t' 

p± 

r^{k'^;z^) = {k'jU^lz^r'U.t" 
PI 

r'yHk'^;z^) = (fcl|*a«[/^|zi)^"C/,i" (D6) 

P± 



RUk'±;z^) = (27i(fcl;z^) + 



{k,k')j 



^Ry,{k'^:z^) + ^Ry,ik'^;z^) V^U^t" (D7) 
ps as J{/ 



2P2 



The corresponding quantities Ryf and ryt are ob- 
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k_/k' 



vertex inside the shock wave. Consider, for example, the 
first diagram in Fig. El Similarly to Sect. I(]3cl we 
calculate the integral over /3' (the p2 component of vector 
k') by taking residues. However, the integral over 
becomes divergent if one takes the term ~ (3' in the three- 
gluon vertex. To deal with such divergence, we must 
retrace one step back and write down the classical field 



^(2 



the form ifHTTIl 



dy{k\ 



p2 



lepo 



l(l)l6 



2D[3A),'>]''{z) (DIO) 




k 



7; 



FIG. 13: Gluon emission by two Wilson lines, 
tained by substitution Uzt'' t^Ul: 



-[k'^\dlu\\z^rh'ul 



{k'^\u—\z^Th'ul 

(fc^|i9«C/^|zi)''™i"C/] (D8) 
PI 



and 



{k,k')^ + ^~tVg 



^Ry,{k'^;z^) + ^Rv2{k'^;z^) y^U^t" (D9) 
(is as ie,/ 



2P2 



Note that Ut"- and t''W carry the independent indices 
of the Wilson lines [oopi,— oopi] and [— oopi,oopi] We 
do not display the color indices of Uzt^ and t''Ul - they 
are always assumed, like (...)([/2t'')J (...)(t^[/])f . Also, 
the formula (ID5II will hold true for Wilson lines in the 



fundamental representation provided one replaces {Ut )j 
by (T^C/)"" and by (r^[/)t™". 

There is a subtle point in the calculation of diagrams 
in Fig. 1131 related to the existence of a term with gluon 



By the equations of motion, one can replace D,Ap 
the r.h.s. of this equation 



(1) 



P.A 



(1) 



2a' s 



T(l) 



2i 



-G,pAi 



(Dll) 



The first term in the r.h.s. of this equation does not 
produce any divergency in P' and can be calculated by 
taking residues. The second term is a contribution with 
the point y (position of the three-gluon vertex) inside the 
shock wave as shown on the last diagram in Fig. H13|l . 
Such terms with the three-gluon vertex inside the shock 
wave are calculated using the formulas for the propagator 
with the initial (or final) points in the shock wave: 



iy \ D2 1^) = [?/*>-oo].y(y|-2 



P^ + iepo ' ' ' "^"'p^+iepQ 

{k\ „2 \y) = [oo,y*],Jfc|^i— |y) (D12) 



p2 + iepo 



p'^ + iepo 



Summarising, j3' in the three-gluon vertex must be re- 
placed by (i — (3 by )^^_Jg^ , and t the difference 
must be taken into account as the term with the gluon 
vertex inside the shock wave. It is worth noting that 
the contribution of the last diagram in Fig. E| (with the 
gluon vertex inside the shock wave) is essential for the 
gauge invariance of the Lipatov vertex (cf. ref. [3^ '). 

The classical field due to the two Wilson lines 

[oopi, — cxjpi]^ [oopi,—oopi\z' is proportional to 

by change of sign 
-> Rv^{ki_ - k'J^, 
'^y. Similarly, 



(2, 

Rvv{k±; z±, z'^_) obtained from (|D5I 
and the replacement Ryt(k — k'Y 
rvtikj_ - k'^Y rvfiki_ - . 



R 



-(X)pi,0Opi]. 



z'j) for the 
[-oopi,oopi]^/ 



classical field due 
is obtained from IID5I) 



to 
by 



change of sign and the replacement Ry{k'^) 

~^ ^yt {k'±)- The Lipatov vertex in the P2^^ — 
gauge takes the form: 
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(L^yl,^{k±;z±,z'^)) = lim k'^{A°;{k)[oopi,-oopi]z [-oopi,oopi]z') 



(D13) 



(fk^ 

"8^ 



(fc- fc')i 



L«''^(fcl;z^)i«f(fc-fc^;zi) 



(fc- fc') 

k[k\ — ki{k,k')_L i^ . . 



kl 



V^'"\k'^-z^)Rly,{k-k'^-z'^) 



kik'iik-k'Y^ 

k. 



-j^R\y{k'^-ZL)R'iv^{k-k'^-z'^) 



-^n2y[k^,z±)+\^gii + Z— 



+ 2r^(fc^;z^) 



-^Rlv^k - fc^; z^) + (^.g^i + 



fci(fc k')i\ I 



kl 



rfy,{k-k'^;z^) 
I 



Again, the Lipatov vertex Lyy^{kj_; z±, z'^) 
of the gluon emission by two Wilson hnes 
[oopi, — oopi]^ [oopi, — oopi]^' is obtained from ljD13|) 
by change of sign and replacement of Ly^{k± — k'^), 
R'y^{kj_ - fcl) and r-y.ik^ - k'^) by L^(fci - fcl) 
Rv{k± — k'j_y and ry{k± — k'^_Y , respectively. The 



vertex Ly\y^^^{k^] z^, z'^) of the gluon emission by two 
lines [—oopi,oopi\z [— oopi, oopi]^/ is obtained from 
(jPlSp by change of sign and replacement of Ly(fc^), 
Ry{k'j^) and ry{k'j_) by the corresponding vertices 
Lt,t(fcl) R'yAk'J and 4^(^-1- " ^1)- 
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Strictly speaking, one should add to the operator the second variational derivative of the source lA9t which contributes 
to the d\U'^U term in the Green function, see ref. jSCjl 

Strictly speaking, to get the effective action Il62t we need only the first term of the expansion of the exponents in Eq. 
1641 1 in powers of A. However, in order to reproduce the full non-linear equation II53II we need the gluon-reggeization terms 
coming from the second order in expansion in A. Formally the gluon reggeization exceeds the accuracy of the semiclassical 
calculation of the effective action; however, when Vi is not large the gluon reggeization is of the same order as 1621 1. 
[37] The requirement of absence of pure gauge fileds at t = oo differs from the condition y4(0, xx) = adopted in the papers 
^ and therefore the fields 1751 differ from ref. For the same reason, the boundary condition ai,2 = a±_ from ref. Q is 
not satisfied by our fields 1751 



